Detecting the Majorana fermion surface state of 3 He-B through spin relaxation 
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The Majorana fermion, which can be useful for topological quantum computation, has eluded 
detection. The 3He-B, recently shown to be a time-reversal invariant topological superfluid, has a 
gap less Majorana fermion surface state. We show here that an electron spin relaxation experiment 
can detect this surface state - its Majorana nature through the Zeeman field direction dependence 
of the relaxation time 1/Ti oc sin 2 where is the angle between the field and the surface normal. 
We propose an experiment setup where an electron inside a nano-bubble is injected below the 3He 
liquid surface. 



Recent development has secured for the Majorana 
fermion a central place in wide range of theoretical 
physics pQ. The chief characteristic of the Majorana 
fermion is that it has only half the degree of freedom as 
the usual complex fermions. It is due to this characteris- 
tic that if neutrions are Majorana fermions we can have 
neutrinoless double beta decay. In recent years, there has 
been great interest in condensed matter systems where 
Majorana fermions can arise. Systematic understand- 
ing of such systems has been obtained through inves- 
tigating their topological properties, which were shown 
to be analogous to those of topological insulator (TI) 
13 HIS EE H E]. Similar to TI, topological supercon- 
ductors (SC) or superfluids have a full pairing gap inside 
the bulk, but have protected gapless state at the edge or 
on the surface. One example is the weak-pairing phase 
of two-dimensional (2D) spinless chiral SC with p + ip 
symmetry [9j [TUJ [TTJ [12] . This system breaks time rever- 
sal symmetry, and can be understood as the SC analogue 
of the quantum Hall (QH) state. The main difference is 
that the chiral edge state of the chiral SC consists of Ma- 
jorana fermions rather than complex fermions as in the 
v = 1 QH state, and thus contains only half the degrees 
of freedom. In addition it was shown that a Majorana 
zero mode is trapped in each vortex core [11] , leading to 
the non-Abelian statistics of vortices [13]. More recently 
the time-reversal invariant (TRI) SC has been proposed 
[14j [15] and classified p~6j H7] . Such topological SC or su- 
perfluid states in two and three dimensions are the anal- 
ogous to the TRI quantum spin Hall (QSH) or the TI 
state discovered recently [I4j [I5j [16]. So far, the only 
definite candidate for the 3D TRI topological SC state 
is the 3 He-B phase [HI [15j [16] , the topological invariant 
of which was first pointed out in Refs. [I2j [18]. In fact, 
the Bogoliubov-de Genne (BdG) Hamiltonian for 3 He-B 
phase is identical to the simplest model Hamiltonian of 
the 3D TI [6j [19], giving rise to a single surface state 
described by the Hamiltonian 



H 



surf 



V F (T ■ (Z X p), 



(1) 



where p is the in-plane momentum, z is the surface nor- 
mal, and <j is the dimensionless spin operator. Despite 
having the same Hamiltonian, the surface state of the 



3 He-B phase consists of a single Majorana cone which 
has only half the degrees of freedom as the surface state 
of the TI which consists of a single Dirac cone. 

There has been recent experimental efforts to detect 
the surface states of 3 He-B [20] [2T] . Despite results con- 
sistent with the existence of the gapless Andreev bound 
state at the surface, these experiments were done on a 
'rough' surface and did not directly detect the Majorana 
cone or the surface state degree of freedom being half 
that of the usual complex fermions. We need a probe for 
a free surface to detect the Majorana nature of the sur- 
face mode, i.e. an analogue of neutrinoless double beta 
decay. There are restrictions on external perturbations 
which can couple to the Majorana surface state of 3 He-B; 
indeed as they are due to the halving of the degrees of 
freedom these restrictions are probably the most distinc- 
tive features of the surface state. The material properties 
of 3 He-B, mainly its very low energy scale, impose further 
constraints on possible experimental methods. Nonethe- 
less, we find that the Majorana nature of the surface 
mode gives rise to some striking and qualitatively dis- 
tinct experimental signatures. 

Surface state of Majorana fermion: First, we show the 
basic similarity and difference between the surface modes 
of 3 He-B and the simplest 3D TRI TI ©Eg. The 3 He-B 
surface mode is derived from the BdG Hamiltonian, 
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where we have used the basis ^BdG( r ) = 
[^^(r), ^^(r), ^l^(r), ^jL(r)] T with the spin quan- 
tization axis along the x-axis (up to rotation by the 
Leggett angle [22j [23] around the surface normal z). 
e p = p 2 /2m is the free fermion Hamiltonian, Ep is 
the 3 He atom Fermi energy, and p± = p y ± ip z . As 
noticed in Ref. [14] [15], H-BdG is formally identical to 
the simplest model of TRI TI with the surface state 
consisting of a single Dirac cone [6] [19]. In both cases, 
the momentum-dependence of the off-diagonal term 
leads to gapless modes bound to the surface [24] [25] . 



The coupling of the spin and orbital degrees of freedom 
for the surface state of the 3D TRI system can be 
understood simply by setting p x = in Eq.Q and 
reduce the system to a 2D TRI system, described by the 
QSH model of Ref. [2 j. This enables us to see that when 
the parallel momentum is aligned along ^/-direction, the 
quasiparticle spin is polarized in ±x direction and the 
— )-spin surface quasiparticle will have dispersion of 
E = —(A/kF)k y (E = (A/kF)k y ). Due to invariance 
with respect to simultaneous spin and orbital rotation 
around z, this coupling of orbital and spin degrees of 
freedom holds for all directions in the x?/-plane. 

Although the BdG Hamiltonian for 3 He-B phase is for- 
mally similar to the model Hamiltonian for the simplest 
TI [6j [19], the fermionic operators that form the bases 
of the two Hamiltonians are quite different. In 3 He-B we 
have particle and hole excitations rather than conduc- 
tion and valence band as in the TI. Since the spin-triplet 
pairing in 3 He-B implies equal spin pairing, we cannot 
distinguish the particle and hole excitation through the 
spin degree of freedom, and thus the annihilation op- 
erator of the negative energy state is equivalent to the 
creation operator of the positive energy state. 

This Major ana nature of the 3 He-B surface mode im- 
poses strong restriction on its interaction with an ex- 
ternal perturbation. To see how this restriction comes 
about, we need to examine the full mode expansion of 
fermion creation and annihilation operators near the sur- 
face. We impose the boundary conditions that the sur- 
face modes vanish at the surface z = and decay expo- 
nentially in the 3 He-B liquid side (where z < 0) of the sur- 
face, albeit much slower than kp- Since the wave vector 
parallel to surface kn remains a good quantum number, 
to satisfy these conditions the surface modes needs to be 
proportional to e* k H" r n sin(k±z)e KZ , where k F = &| + k\ 
and k > 0. Inserting this to the BdG equation Eq.([2| 
gives k = A/hvp and reduces Eq.Q to an effective sur- 
face Hamiltonian of Eq.Q. Therefore, for our surface 
mode expansion we use the result from the TI but also 
take into account the artificial doubling mentioned above: 
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where </>k = a.Tctaxi(k y / k x ) and is a normalization con- 
stant of the mode k (see SOM for details). Note that once 
we ignore the gapped modes (eigenenergy greater than 
A), we obtain the Majorana condition ^_>(r) = ^l^(r) 
and V^(r) = ^jl(r). What this means is that the 
local creation and annihilation operators for a fermion 
with its spin polarized parallel to the surface is in- 
distinguishable once we ignore modes with eigenenergy 



greater than A, thus reducing the degrees of freedom 
by half. Instead of the usual fermion anticommuta- 
tion relation, these Majorana operators would form Clif- 
ford algebra, ^2 a {^a( r )^ ( r/ )} = 25(r — r') (where 
<j, a' =— <—). It follows that it is impossible to construct 
the spin-polarized local density p a (r) = ^(r)^ a {r) out 
of the gapless modes if the polarization axis is parallel 
to the surface. This means that with the gapless sur- 
face mode, we can neither construct the local density 
operator p(r) = ^2 a ^i(^)^a(' r ) nor the components of 
the local spin density operator parallel to the surface, 
I x = (^^_-0t_^)/2 and 4 = ($L^<- + $1^-0/2. 
However, it is possible to construct the component of spin 
density operator perpendicular to the surface, I z (r) = 
— z , 0_>(r)'0<_(r) [10]. So in 3 He-B the surface state does 
not contribute to the local density fluctuation while its 
local spin density is effectively Ising for T <C A, which 
means that the local external perturbation can excite the 
surface state only if it couples to I z ; this is a direct con- 
sequence of the halving of the degrees of freedom. 

Therefore, to detect the surface state and its Majo- 
rana nature, it is best to measure dynamic susceptibility 
arising out of these gapless modes. From the discussion 
above we see that the dynamic spin susceptibility tensor 
of the surface state has only single nonzero component: 
X zz , which we can calculate from Eq.(|3|. Anisotropy 
this drastic cannot be obtained from spin-orbit coupling 
of the complex fermions such as we see in the TI sur- 
face state. So we conclude that the resonant spin spec- 
troscopy is the best probe for the Majorana surface mode. 
The extreme anisotropy of the spin susceptibility should 
be revealed through striking anisotropy in the spin spec- 
troscopy. Due to the gapless dispersion, there will be no 
g— a/t SU pp ress ion of this anisotropy. We now need a 
spin probe that best fits the material property of 3 He-B. 

ESR - spin spectroscopy: We propose electron spin re- 
laxation (ESR) as the best spin spectroscopy on the 3 He- 
B surface state. Our basic idea is to introduce some ex- 
tra electrons to 3 He-B, apply a weak DC magnetic field 
(which satisfies H <C T//i B ; note A/fi B « 26.2G [26]), 
excite the electron spins through resonance, and then let 
these electron spins relax through interaction with the 
surface state. This relaxation process would probe the 
dynamic spin susceptibility of the 3 He-B in a way analo- 
gous to the way the nuclear magnetic relaxation (NMR) 
is used to probe the dynamic spin susceptibility of elec- 
tron in a crystalline system. Such probe should reveal 
the drastic anisotropy of the dynamic spin susceptibility 
of the surface state due to its Majorana nature. More 
explicitly, we start from the spin relaxation rate formula: 
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where P(q,z e ) is the static form factor of the electron 
(obtained from Fourier transforming the xy coordinates 
of the probability density of a single electron), A+ is the 
component of the interaction that flips the electron spin 
with respect to the direction of the Zeeman field, z(z') 
and z e (z' e ) are the ^-coordinates of the 3 He atoms and 
the electron respectively, and uol = gfis/h is the Larmor 
frequency of the electron. This formula would look like 
the standard NMR relaxation formula [27] if we drop out 
the z dependence, the electron form factor P, and restore 
the isotropy of the dynamic spin susceptibility. Eq.Q 
implies the dependence of 1/Ti on the direction of the 
Zeeman field, because A+ couple I z to the component of 
the electron spin perpendicular to the Zeeman field. 

To illustrate this dependence on the Zeeman field direc- 
tion, we consider a simple contact interaction model for 
the coupling between the electron and 3 He atom spins. If 
we set the magnetic field direction as z / = z cos 6+ic sin 
we can write down the contact interaction as H contact = 

-A con tactIzS z = ~ A contact I z [S z > COS 0- \ (S++SL) SHI 0] , 

giving us A + = A contact sin#. Inserting this into Eq.Q, 
we obtain 1/Ti oc sin 2 0. In other words, the electron spin 
does not relax at all for perpendicular field! By contrast, 
the same model gives us 1/Ti independent of for the 
surface state of the simplest TI, q summation canceling 
out the spin susceptibility anisotropy. 

Realistic calculation can still give us this drastic 
anisotropy of spin relaxation. In 3 He-B, the main chan- 
nel of spin-spin coupling is the dipole-dipole interaction, 
mainly because an electron strongly avoids contact with 
3 He atoms. With the dipole-dipole interaction, we do 
have coupling between I z and S XjV : 
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where 7 is the gyromagnetic ratio of a 3 He atom and g is 
the Lande ^-factor of an electron. However, for the elec- 
tron below the liquid surface, the S X:V terms of Eq.([5| 
may have little effect; because z > for helium atoms 
'below' the electron and z < for helium atoms 'above' 
the electrons, the coupling to S XiV from the helium atoms 
above cancels out the coupling to S XiV from the helium 
atoms below. Since the spin interaction is effectively Ising 
(that is, Hd oc —I z S z ), we have 1/Ti oc sin 2 0, as we ar- 
gued the previous paragraph. By multiplying sin to the 
2D Fourier transform on the coefficient of the I Z S Z term 
of Eq.(|5]), we obtain A+faz) = - ^9^h qe - q \ z \ gin6 > 
As the next step, we need to devise an experimental setup 
to relax the electron spin by the 3 He-B surface state. 

Electron bubble: A crucial constraint on the relaxation 
rate is how well the electron is localized. Whereas in 
the NMR, we can assume that a nucleus is a point-like 
object, we cannot make the same assumption for elec- 




FIG. 1: Illustration of the surface state of 3 He-B phase con- 
sisting of a single Majorana cone, where the E < part of 
the quasiparticle spectrum is redundant and indicated with 
the dashed boundary. Also shown are the dimensions of the 
bubble electron when we apply a perpendicular electric field 
of 150V/cm. Note how small the size and depth of the bubble 
is compared to the depth £ of the surface state, for which we 
take the weak coupling approximation Hvf/^ as in Eq.Q. 

trons in ESR and hence the introduction of the static 
form factor P(q) in Eq.Q. Due to the Heisenberg un- 
certainty principle, the more delocalized the electron is 
in the real space, the more rapidly P(q) falls off with 
q. This suppresses the spin relaxation for processes that 
result in a large momentum change for 3 He atoms and 
hence suppresses 1/Ti. For this reason, 1/Ti is very small 
for an electron sitting on top of the 3 He liquid surface. 
Even when electrons above the surface form a Wigner 
crystal, the zero-point displacement is greater than 10% 
of the lattice constant for the lattice constant < 1/im 
[28] . There is a limit to reducing the lattice constant 
as we need to keep the dipole-dipole interaction between 
adjacent electrons much weaker than the interaction be- 
tween an electron and 3 He atoms. In order to enhance 
the electron localization significantly, we need to place 
the electron under the 3 He liquid surface. 

Once it is injected below the 3 He liquid, an electron 
settles into a well-localized metastable state below the 
surface. It cannot be easily ejected from the liquid due 
to an electrostatic energy barrier at the surface arising 
from the fact that helium is dielectric [29]. From the 
electric field boundary condition E z \ z= o+ = €dE z \ z= o- 
we see that the polarization of atoms causes the surface 
of dielectric to repel an electron below the surface. By 
tuning an electric field perpendicular to the surface, we 
can adjust the equilibrium distance \b\ between the elec- 
tron and the liquid surface from right near the surface, 
~10nm, to below the surface state, \b\ > £ [29] [30]. Be- 
low the liquid surface, an electron opens up a nano-sized 
cavity and becomes trapped inside of it to avoid the en- 
ergy cost due to the negative electron affinity of helium 
atoms. The size of this 'bubble' is determined by compe- 
tition between the zero-point kinetic energy of the con- 
fined electron Ezp = h 2 / (8mR 2 ) and the surface energy 
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FIG. 2: The electron spin relaxation rate (in 10 _3 Hz) due 
to the surface state through dipole-dipole interaction for the 
magnetic field applied parallel to the surface. From the top to 
the bottom curve, the applied perpendicular electric field is 
150V/cm, lOV/cm, and 1.5V/cm respectively, which gives us 
the bubble depth |6|= 22.5nm, 87.4nm, and 225. 2nm respec- 
tively (for consistency with Eq.([3]) the depth of the surface 
state is set £ = /i^i? / A 237nm ). 

of the cavity Es = 47ri? 2 a, where R is the cavity radius 
and a =0.156 erg/cm 2 is the surface tension of the he- 
lium liquid [31]. This gives us the electron localization 
R = [h 2 / (327rma)] * =2.35nm, far better than what we 
obtain above the surface. Fig. 1 shows this electron bub- 
ble radius compared to the depth of the surface state. 

Our ESR rate calculation shows signatures of both the 
Majorana nature and gapless dispersion. For electron 
bubbles placed at 22.5nm, 87.4nm, and 225. 2nm below 
the surface, we find that the relaxation rate is 5 x 10 2 
times faster for the parallel field (9 = tt/2) than for the 
perpendicular field (0 = 0), implying that we effectively 
have 1/Ti oc sin 2 relation. As shown in Fig. 2, for 
the bubble depth of 22.5nm, the relaxation rate 1/Ti is 
approximately 10 3 sec (see SOM for details) The absence 
of the exp[— A/T] suppression in 1/Ti versus T behavior 
characteristic of the bulk quasiparticle is the consequence 
of the gapless dispersion on the surface. However the 
relaxation rate anisotropy will be reduced if we include 
contribution from bulk condensate, which has isotropic 
nonzero spin susceptibility [22j [23] . 

In conclusion we have proposed a realistic experiment 
setup to observe the Majorana fermion surface states of 
the topological superfluid 3 He-B phase. Due to the Ma- 
jorana nature of the surface state, the spin density op- 
erator is purely Ising-like, polarized perpendicular to the 
surface. Through an ESR experiment, we can show both 
gapless dispersion and extreme anisotropy of the dynamic 
spin susceptibility. Our experimental setups for the ESR 
measurement uses electron nano-bubbles placed below 
the liquid helium surface, giving rise to the 1/Ti oc sin 2 
dependence on the magnetic field direction. Such a di- 
rect experimental observation of the Majorana fermion 
would enhance our fundamental understanding this ex- 
otic particle and the nature of the topological superfluid, 
and pave the way for topological quantum computing. 
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Notes added: Near the completion of this work, we 
learned that Nagato et al. independently obtained the 
spin susceptibility anisotropy of the surface state [32] . 
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THE SUPPORTING ONLINE MATERIAL 

Surface state mode expansion 

In solving the BdG Hamiltonian Eq.(2), we make the 
following weak-pairing approximation: 
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Now when the in-plane momentum is aligned along y- 
direction, other than the quasiparticle spin-polarization 
in ±x' direction, it is essentially identical to the 2D SC 
with Py+ip z pairing for | —>—>) pairs andp y — ip z pairing 
for | <— <— ) pairs. From the edge state of such SC, we can 
see that with the in-plane momentum aligned along y- 
direction, — )-spin surface quasiparticle will have dis- 
persion of E = — (A/k F )k y [E = (A/k F )k y ). Such align- 
ment of orbital and spin degrees of freedom should hold 
for all direction in the x?/-plane. Therefore the full mode 
expansion gives us Eq.(3). Note that neither introducing 
anisotropy between the in-plane (p±A/p F ) and perpen- 
dicular (p z A/p F ) components of the gap nor taking into 
account the possible z dependence of the perpendicular 
component of the gap is going to change the mode ex- 
pansion qualitatively. 



anomalous part of the Green function do not contribute. 
This means that \ zz is the same as the 3D strong topo- 
logical insulator surface state with the chemical potential 
at the Dirac point. Therefore, the imaginary part of the 
dynamic spin susceptibility in this limit is 
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Lastly, there is the static form factor of the electron 
P(q,z). This is a 2D Fourier transform of the modulus 
square of the single electron ground state wave function. 
In the case of the electron in a bubble, we can approxi- 
mate the electron wave function to vanish at the bubble 
boundary, so we can set for the modulus square of the 
ground state wave function 
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Calculating the relaxation rate 

To calculate the spin relaxation rate Eq.(4) due to 
the dipole-dipole interaction of Eq.(5), we need to cal- 
culate from Hjj the matrix element that couples to 
S- = (S x cos — S z sin 0) — iS y in the momentum space. 
To obtain this, we define A^(r) from Hp = I z (t)Ad(t — 
r e ) • S(r e ). 2D Fourier transform gives us 
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From this, we can see that the component of 
couples to S- is 
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where R is the bubble radius, if we take the center of 
the bubble to be the origin. A 3D Fourier transform of 
P(r) approximates to exp[— (R 2 /2it 2 )(q 2 
can make approximation 
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In calculating the imaginary part of the dynamic sus- 
ceptibility, we take the uj f — > limit. In this limit, the 



where we now take the coordinate of the center of the 
bubble to be (0, 0, b). We also note that since the energy 
difference between the electronic ground state to the first 
excited states is much larger than the pairing gap of 3 He- 
B. 

With a given bubble depth |6|, we can calculate the 
ESR relaxation rate from Eqs.([8|, and (JTTJ . We find 
• that the relaxation rate due to the helium atoms in the 
(3) \z — b\ < R region is only about 4% of the contribution 
from the \z — b\ > R region. Therefore we only give the 
result for the latter region. The result can be expressed 
in most part by dimensionless integrals. For the Zeeman 
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field parallel to the surface (6 = n/2) we obtain 
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For helium-3 liquid e d = 1.04 [fj; this gives \b equi i\ 
=22.5nm for E = 150V/cm. We have simply set b = 
b equ ii and have not considered any fluctuation around 
this equilibrium point. 
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(12) 



where x = k/kp, V = g/2fc and T = ksT/A. For the 
perpendicular (0 = 0) Zeeman field, we obtain the relax- 
ation rate 
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ex/T)2j Q dy 



[erf (tt/V2 - V2xyk F R/n) 
erf (tt/V^ + V2xyk F R/ir)] 2 



2R 

e € 



e -2xyk F (b-R) 



_ 2fi 



xyk F ^ - 1 
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We see in Eqs.([12| and ([13]) that sgn(z) in Eq.([8]) leads 
to the cancelation between contributions from —b + R< 
z < and from z < —b — R. For fci? = 7.88/nm, £ = 
237nm, i? = 2.35nm, Eq.Q is larger than Eq.Q by 
thee orders of magnitude at \b\ = 22.5nm, 87.4nm, and 
225. 2nm; the dependence on \b\ is quite weak. Fig. 2 is 
the plotting of Eq. (12) for different values of \b\. 

Finally, we note that the relation between the applied 
perpendicular electric field and the bubble depth is 



equil | 



AttsqE Ae d {e d + 1) 



1/2 



(14) 



This comes from noting that a charge q below dielectric 
(ratio Ed) surface at this depth (z = 6; note b < 0) in- 
duce charge on dielectric surface through polarization. 
This surface charge is effectively equivalent to having an 
image charge q(sd — l)/(£d + 1) at z = —6, which means 
that the dielectric surface repels the charge below the 
surface. Therefore, when we apply the constant perpen- 
dicular electric field E, the net force the electron bubble 
feels is 
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F z {b) = 



4^o(26) 2 £ d (^ + l) 



eE. 



(15) 



6 



